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ABSTRACT 



We generalize the derivation of dynamo coefficient a of Field et al. (1999) to 
include the following two aspects: ffist, the de-correlation times of velocity field 
and magnetic field are different; second, the magnetic Prandtl number can be 
arbitrary. We find that the contributions of velocity field and magnetic field to 
the a effect are not equal, but affected by their different statistical properties. 
In the limit of large kinetic Reynolds number and large magnetic Reynolds 
number, a-coefficient may not be small if the de-correlation times of velocity 
field and magnetic field are shorter than the eddy turn-over time of the MHD 
turbulence. We also show that under certain circumstances, for example if the 
kinetic helicity and current helicity are comparable, a depends insensitively 
on magnetic Prandtl number, while if either the kinetic helicity or the current 
helicity is dominated by the other one, a different magnetic Prandtl number will 
significantly change the dynamo a effect. 

Subject headings: Dynamo, MHD turbulence 
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1. Introduction 

Dynamo theory is an attempt to understand tlie process of magnetic field generation by 
self-inductive action in electrically conducting fluids. The theory has been used to explain 
the magnetic field generation in many celestial objects such as the Sun, the Earth and the 
Galaxy. Dynamo theory divides into a kinematic regime and a dynamic regime. In the 
kinematic regime, the velocity field is prescribed, and dynamo theory studies the physics of 
magnetic field under the determined velocity field. In general, kinematic dynamo theory 
focuses solely on the mathematical solution of the induction equation, 

— = Vx(VxB) + AV2B, (1) 

where V, B are the velocity and magnetic field, respectively, and A is the magnetic 
diffusivity. 

In the dynamic regime, the velocity field can be modified by the magnetic field, so 
that the dynamic dynamo theory must then consider both the induction equation and the 
momentum equation 

^ = -V ■ VV + B ■ VB + z/V^V - VP + f , (2) 

ot 

where u is the molecular viscosity, P the total pressure, and f the external forcing term. 
The mean-field electrodynamics(MFE) developed by Steenbeck, Krause and Radler(1966) 
through the two-scale approach provides essential insights into the relation between the 
statistical properties of turbulence and dynamo effects, namely, the a and j3 parameters 
of dynamo theory. Let the large-scale magnetic field and velocity field be B and V, and 
the fluctuating components of magnetic field and velocity field be b, v, respectively. The 
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two-scale separation of MFE gives the equation for B asQ 

dB 



dt 



V X (V X b) + Vx < vx b > +AV^B. (3) 



The so-called turbulent electromotive force, £^=<vxb>, is related to the dynamo a and 
P effects through 

S = aB + (3VxB. (4) 

MFE in the kinematic regime gives (see section 7.3 of Biskamp, 1993, and references 
therein) 

a = -- < V- Vx v>, /3 = - < v^ >, (5) 

where r is the velocity de- correlation time. Largely as a result of the development of MFE, 
the kinematic aspect of dynamo theory has been broadly understood. Several monographs 
have been devoted to this subject (Moffat 1978, Krause and Radler 1980, Zeldovich et al. 
1983). 

But the nature of dynamo theory in the dynamic regime is still in debate. The back 
reaction of the magnetic field on the velocity field will modify the expressions for the 
dynamo a and (3 coefficients. The numerical simulation by Pouquet et al. (1976) was among 
the first to point out that with back reaction, to lowest order in B the dynamo a effect 
should be modified to 

ttfer = a^ + ai, (6) 

where a„ is the a in (5) and a^ is proportional to < b ■ V x b >0. This new term will 
reduce the classical kinematic dynamo a effect to a certain degree. Some early criticism 



^< ■ > and ~^ are interchangeably used throughout this paper to denote ensemble average. 
"^abr was called residual torsality by Pouquet et al. Following Field et al. (1999), we call 
< b ■ V X b > current helicity in this paper. 



-5- 

of MFE was discussed by Piddington (1970, 1972abc, 1975ab). He argued that kinematic 
solar dynamo theories do not account for the removal of the large amounts of flux generated 
each solar cycle. Recent objections to dynamo action have their root in the problem of 
small-scale magnetic fields. For astrophysical systems such as the Galaxy, where magnetic 
Reynolds number Rm = Lvq/X is large, some authors(Cattaneo and Vainshtein 1991, 
Vainshtein and Cattaneo 1992, Gruzinov and Diamond 1996) argue that the magnetic 
energy at small scales, < b^ >, is much greater than the magnetic energy at large scales. 



2 

B , through the relation 



< b^ >= R^b\ (7) 



According to Cattaneo and Vainshtein (1991), the magnet ohydro dynamic turbulent dynamo 
will stop operating as soon as relation (7) is obtained, a process that can happen in much 
shorter time than turbulent eddy turn-over time (Kulsrud and Anderson, 1992). Gruzinov 
and Diamond (1994) base their argument on the conservation of squared vector potential 
for 2D MHD and magnetic helicity for 3D MHD, and claim that dynamo a effect will be 
quenched in systems of high magnetic Reynolds number as follows 



a 

a = 



1 + RjB/vof 



where a^ = — t/3 < v ■ V x v > is the classical result in (5). The numerical simulation 
with periodic boundary conditions by Cattaneo and Hughes (1996) supports relation (8) for 
the particular value Rm = 100 and various values of B. Relation (8) is completely different 
from a previous estimate of a effect made by Kraichnan (1979). He argued that even in the 
high magnetic Reynolds number limit, the a effect will not be quenched. Rather, it has the 
following relation with the classical estimate, a^, and B, 

" ;= — 72- (9) 



l + {B/voy 
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Kulsrud (1999) questioned the derivations of Gruzinov and Diamond (1994, 1995, 
1996) by arguing that one of the their resuhs, a = /9B ■ V x B, for large Rm, leaves out 
the contribution from v ■ V x v completely. Note that Gruzinov and Diamond derived (8) 
from the conservation of magnetic helicity, A ■ B, which is based solely on the induction 
equation. In the dynamic regime, a velocity field that is constantly driven by external 
force will be modified by the growing magnetic field, a process that cannot be wholly 
understood by using only the conservation of magnetic helicity. Leaving this process out of 
the discussion of dynamo action is questionable, and a more complete account of dynamo 
theory in the dynamic regime must also consider momentum equation and the role of 
external forcing terms in maintaining the turbulence. Field, Blackman and Chou (1999) 
considered a simplified model of MHD turbulence. The external forcing term in their model 
has the freedom to drive a MHD turbulence that is independent of the presence of any 
large-scale magnetic field. The turbulent velocity field and magnetic field were treated on 
an equal footing. Their result on the a dynamo effect depends on the statistical properties 
of MHD turbulence that are independent of B, and agrees approximately with that of (9) 
but disagrees with (8). 

In this paper, we generalize the work by Field, Blackman and Chou (1999) in the 
following two aspects: First, we relax the assumption made by Field, Blackman and Chou 
that the de-correlation times of velocity field and magnetic field are the same. Second, we 
relax the assumption in Field et al. that magnetic Prandtl number, Pr = u/X, is one. We 
also give brief discussions on how the nonlinear interaction between B-induced quantities 
(v' and b' in Field et al. ) changes the final result for a, and compare our result with (8) 
and discuss where the difference originates. This work can be considered as a generalization 
of Field, Blackman and Chou (1999), yet itself is complete and self-contained. To simplify 
matters, we discuss only the dynamo a effect (not (3 effect) in this work by assuming a 
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constant large-scale magnetic field, B. 

The structure of this paper is as follows: in section 2, we discuss our model based 
on the work of Field et al. In section 3 and appendix A, we give our derivation of the 
coefficient of dynamo a effect. In section 4, we talk about the dependence of the a effect 
on the kinetic Reynolds number Rf,, the magnetic Reynolds number R^ and the magnetic 
Prandtl number Pr{= v/X) based on our derivation in section 4. In section 5, we discuss 
possible modifications to our result in the presence of strong nonlinear interaction between 
B-dependent components of the MHD turbulence, and compare our result with previous 
derivations of the a effect. Conclusions are made in the final section, section 6. 

In Table 1, we list the notations used throughout this paper. The physical meanings of 
these quantities are also explained briefiy in this table. 

2. Separation of v', b' from v*^*^) and b*^°^ 

As in the work of Field, Blackman and Chou(1999), we consider only incompressible 
fiuids. In our model, we first distinguish four different scales in the system. We denote 
the size of turbulence energy containing eddy as /q, and denote the dissipation scale of 
the turbulence as Id- For turbulence of large kinetic Reynolds number, we have /q ^ Id- 
An "ensemble average scale", denoted by L ^ /qj is used to carry out the calculations of 
averaged quantities, < ■ >. Finally, we denote the scale of the whole physical system as 
S*, which is the typical scale for the variations of averaged quantities < • >. So we have 
the relation S^L^Iq^Io- In our model, there are two large-scale quantities, B and 
V. We assume both of these quantities are constant, i.e., 5*^00. We set our reference 
frame to that moving at V and henceforth omit terms of V. The presence of non-zero B 
is related to the a effect in the following two aspects: first, B can be amplified due to the 



dynamo a effect; second, for large B, a effect will be quenched. 

Next, we make the distinction between the B-independent components, v*^*^) and b*^°) 
from the B-induced components, v' and b'. In the absence of a large-scale magnetic field 
B, we assume that a homogeneous, isotropic, steady-state MHD turbulence be maintained 
by an external force, which operates at scale ~ /q- To simplify matters, we consider only 
the case that the external force is homogeneous, isotropic and helical. We further assume 
that within the scale range, [e/o,e^^^o] where e > 1, the turbulence is largely hydrodynamic 
in character. This assumption can be justified if the external force, which is independent of 
the velocity and the magnetic fields, is the only external energy source to drive the velocity 
field. A homogeneous, isotropic, helical external force will constantly drive a homogeneous, 
isotropic, helical velocity field near the forcing scale. The equipartition between the velocity 
field and magnetic field is not achieved within range [e/o,e^^/o]- Rather, the kinetic energy 
within this range surpasses the magnetic energy, a phenomenon that has been found in 
many numerical simulations (see below). From the scale e^^/o down to viscous cut-off scale 
~ Id of the velocity field, equipartition can be maintained according to the MHD turbulence 
model of Kraichnan(1965) or Goldreich and Sridhar(1995). The simulation by Pouquet 
et al. (1976) with simplified DIA equations shows that e ~ 3. Simulations by Kida et al. 
(1991) with a spectral method gives e ~ 5. Both authors obtained homogeneous, isotropic 
MHD turbulence with no large-scale magnetic field present, i.e., B = 0. We denote the 
velocity and magnetic field in this kind of homogeneous, isotropic turbulence as v*^°^ and 
h^^\ and clearly they satisfy the equation 

5,b(°) = V X (vW X bW) + AV^b^o), (10) 

and 

aiv(°) = -v(°) ■ Vv(°) + b(°) ■ Vb(°) + z/VVW - VP(°) + f, (11) 
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where A and z/ are magnetic resistivity and molecular viscosity, respectively, f is the 
external forcing term. 

As we mentioned before, B can be amplified by a. effect, and a growing B will also 
attempt to reduce a dynamo effect. To tackle this problem, we impose a large-scale 
magnetic field, B, on the MHD turbulence of v*^"^ and b^°\ and then follow the dynamics of 
the MHD turbulence system, v'^''^ will stretch B, so that a new component, b', of magnetic 
field that is dependent on B, is generated. The new magnetic field, B + b' + b^^-*, will exert 
Lorentz force on the velocity field so that a new component of the velocity field, v', which 
depends on B, is in turn generated, v^^^ and b*^") still satisfy equations (10) and (11) even 
in the presence of B. This is because the change in the total velocity field, v' = v — n'^^\ is 
due to B, while the component v*^*^) is driven primarily by external forcing term f , which is 
independent of B. Bearing this in mind, and assuming that we have v^°) and b'^'^^ already 
at our disposal^, we write down the relation between the known quantities, {v'^'^^ b'^^-', B}, 
and the unknown quantities, {v', b'}, as 



9iv' = B ■ V (b' + b(°)) + z/VV, (12) 

9tb' = B ■ V (v' + v(°)) + AV'b'. (13) 



(12) and (13) show that the generation of B-dependent component of magnetic field is due 
to the stretching of B by both the newly generated, B-dependent v', and the B-independent 
velocity component v'-''-'. The generation of B-dependent component of velocity field is 
due to the Lorentz force that is exerted by the combination of b' and b*^"-' . In this sense, 
the terms involving v*^'') and b*^") in (12) and (13) can be considered as sources for v' and 



■^By, say, numerical simulations of equations (10) and (11), or observations of MHD 
turbulence with no large-scale magnetic field present. 
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b', thus the nonhnearity of v' and b' will largely come from the nonlinearity of the MHD 
turbulence, namely, v*^'^-' and b^^-*. Therefore, we expect that the statistical properties of v' 
and b' depend on the statistical properties of v'-^-' and b^^-* . 

In deriving equations (12) and (13), we made an important assumption. The validity of 
this assumption has been discussed in Field, Blackman and Chou (1999). The assumption 
is that the ratio of the decorrelation time Tdcor of both magnetic and velocity fields, to the 
eddy turn-over time, Teddy i is small: 

^^ < 1. (14) 

'eddy 

With such assumption, we dropped the following nonlinear terms: b' • Vb', b' ■ Vb(°\ 
b' ■ Vv', b' ■ Vv(o), v' ■ Vb', v' • Vb(o), v' • Vv', v' • Vv^, b^o) ■ Vb', b^o) ■ Vv', v^") ■ Vv', and 
v'^°^ ■ Vb'. These terms are all smaller than the time derivative terms, dth' and dt'v' ■ In fact, 
the ratio of these two terms to the nonlinear terms that we dropped is approximately l£^^. 
We make this assumption because (for more details, see also Field et al. , 1999) numerical 
simulations by several authors (Pouquet et al. 1976, Kida et al. 1991, Brandenberg 2000) 
have shown that the MHD turbulence within range [e/o, e^^^o] for e > 3, is dominated 
by velocity field. In other words, the MHD turbulence is largely hydrodynamic within 
this range, where most of the energy is concentrated. For pure hydrodynamic turbulence, 
experiments (Pope 1994) have shown thatQ 

^^ = 0.2-0.3. (15) 

"Teddy 



^Our recent numerical simulation (Chou and Field, 2000) of 3D MHD turbulence also 
supports the validity of assumption (14). We found that for 3D MHD turbulence in steady 
state (50 < Rm < 300, B < Vrms), Tdcor /Teddy ~ 0.22 for vclocity field and ~ 0.25 for 
magnetic field. Similar results are also obtained by Brandenberg(private communication). 
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Therefore the time derivative terms in (12) and (13) are greater than the nonhnear terms 
that we dropped to get (12) and (13). Note that the assumption of (14) does not imply 
that the turbulence is of only one scale Iq- In fact, the turbulence in this work can be of 
multiple scales. However, if the turbulence at scales < e~^/o is in approximate equipartition, 
its contribution to the electromotive force, < v x b >, will be small because of the Alfvenic 
wave-like motions of the MHD turbulence in the equipartition scales. In other words, 
we believe that the dynamo a— effect to be discussed in the next few sections is mostly 
attributed to the turbulence within scales [e/o, e^^^o] for e > 3, and various numerical 
simulations mentioned above support such picture. 

In this work, the ratio of B to the rms of turbulent velocity, Vrms, can be of any value; 
therefore we keep all terms that involve B in the MHD equations for v' and b'. Indeed, by 
keeping all terms of any order in B/vrms, we include the nonlinearity from the back reaction 
of B to the turbulence. Also, under condition (14), our treatment of a effect does not 
require any constraint on the ratio of B to the rms of magnetic field, brms- Thus the back 
reaction of b on the turbulence is also taken into account in our model. This is different 
from classical dynamo theory that assumes that brms is smaller than B. 

In the following, we assume that equations (10) and (11) have been solved so we have 
v*^*^-* and h^^' at our disposal. In fact, we require only the statistical properties of these 
zeroth order velocity and magnetic fields. Now, equipped with equations (12) and (13), and 
the known quantities v*^°) and h^^\ we will be able to calculate the relation between {v', b'} 
and {v*^''\ b*^°^} and, furthermore, calculate the electromotive force, < v x b > and relate it 
to dynamo a-effect. The calculation is performed in next section and appendix A. 
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3. The Derivation of Dynamo a-coefRcient In the Presence of Nonlinear 

Effects of v(o) and b^o) 

Throughout this paper, we define the Fourier transform F of a function F as 

F(k,cu)=f^) f r^^ F{^,t)e-''^^''-'^'^dtdx; (16) 

\2tt J JvJ~oo 

i.e., we transform F(x, t) both in time and in space. This method is similar to the one used 
by Roberts and Soward(1975), and we assume B = constant throughout this work. The 
equations of b' and v', (12) and (13), are transformed into Fourier space as 

[-lu + ue) v'j = i{k ■ B)(6; + bf), (17) 

[-ZU + \e) hr = i{k ■ B){v'j + vf). (18) 

Here j = 1,2,3. We will use the Einstein summation convention throughout this paper. 
Solving (17) and (18) gives the following relations: 

6;.=a(k,c.)«f -C(k,c.)6f, (19) 

v'^=UKu^)bf-aKu^)uf\ (20) 

where 



i(k-B) 

UK UJ) = '^^^g. — , (21) 

t(k-B) 



1 



and 



UKu^)= -g fe , (22) 

^ + {-iu)+Xk^){-iuj+iyk2) 



C(k, uj) - ^ ^ ^^-- . {26} 



{~iuj+Xk^){-iuj+uk'^) 

The solutions (19) and (20) relate the induced velocity field and magnetic field to the 
background b*^*^) and v*^°^. Except for the condition (14), we set no constraint on these 
zeroth order components, i.e., the components that are independent of B. 
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To calculate dynamo a-effect, we use the following relation 



(v X b)„ = aBn. 



(24) 



Here v = v'^''^ + v', and b = h^^^ + b'. We then apply the properties of the Fourier transform 
and write the correlation between v and b as 



{vi{x,t)bj{x,t)) 



rfkrfo^e^e^-^-^*) 



dk'du' {vi{k', Uj')bj{k -k',uj- uj'] 



(25) 



Let gij = (vi(k',uj')bj(k — k',ti; — on')); then with (19) and (20), we have 



gij = [l-C(k-k',cu-cu')-C(k',cu') 

+ C(k', cu')C(k -k',LU- LU')] P°°(k -k',iu- iu')5ik, iu) 

+ [6(k -k',uj- Lo') - C(k', cu')6(k -k',uj- Lo')] g°°(k -k',uj- uj')6{k, uo 

+ [e.(k', uo') - U^, oo')ak - k', cu - oo')] ^°0(k - k', cu - oo')5{k, u) 

+ [6(k', uj')U^ - k', .; - uj')] P°°(k', uj')6{k, u). 



(26) 



Here we've used the homogeneity and isotropy of v*^°^ and h^^\ and the following definitions 
of correlation functions and their corresponding Fourier transforms: 



t;f°)(x,t)t;f(x',t'))=Q?,°(x'-x,t'-t), 



,(o)r^ +u(o)r^' +'^ 



t;r(x,t)6;.'^^(x',tO)=P°°(x'-x,t'-t), 
6!°)(x,t)6f(x',t')) = 5g°(x'-x,t'-t), 



^(o)^ 



^(o)^ 



1)}"^ (k, uj)v'l'> (k', cu') ) = Q^"(k', cu')5(k + k', cu + cu' 



{);°)(k, Lu)bf\k', iu')) = P°°(k', cu')5(k + k',iu + iu'), 
b'i'\k, u;)bf\k', u')) = Sf^{k', u;')6{k + k', cu + u'). 



(27) 

(28) 
(29) 
(30) 
(31) 
(32) 



-14- 

With such definitions of correlation functions and tlie properties of Fourier transform, (24) 
can be further calculated as (see Appendix A for details) 

(v X b)„ = M„ + iV„, (33) 

where 

M„ = J J dkduj [6(k, u) - C(-k, -uj)UK uj)] e„«,Q?°(k, iv), (34) 

and 

N^ = JJdkdu [U-K -UJ) - ak,uj)U-K -^)] eni^S^'iKu). (35) 

The above integrals can be calculated with the following assumptions about the MHD 
turbulence. Suppose that we can separate variables in the Fourier transformed correlation 
functions and write Qy (k, cj) = X[^°(k)T(c<j) and Sij^{k,uj) = 0^j'(k)J{uj). In this paper we 
also assume that T(ip) = e~l'^l'^'= and J(v5) = e~''^'^'= where Lp = t — t\ uJk and Xk are positive 
and functions of the wave vector, and ujk = l/Tt,_cor and Xk = ^/Tb_cor where Tv_cor is the 
de-correlation time of velocity field at k, and Tii_cor is the de-correlation time of magnetic 
field at k. Again, because of the homogeneity and isotropy, we have that (Krause and 
Radler, 1980, p. 75) 

XfP(k) = A{k)Si, + B{k)hk, + teki,kkC{k), (36) 



and 



Oj;"(k) = D{k)5i, + G{k)kik, + teki,kkH{k), (37) 

where A, B, C, D, G, H are functions of A; = |k|. Here we use C{k) and H{k) to denote the 
helical parts of the spectra of the velocity field and magnetic field, respectively. With these 
assumptions and tedious evaluations (see Appendix A for details), we have much simpler 
integral forms of M and A^ as 



1 r+oo 

M„ = -- / dk 
I Jo 



-^A,{k)U,{k) 



Bn, (3^ 
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where 



iV„ 



1 

+ 9 / 

2 JO 



dk 



2z/ 
Z/+ A 



Mk)Mk) 



Br 



Ab{k) = 



87ik*C{k) 

8nk^H{k) 
Xk 



(39) 

(40) 
(41) 



n,(A;) 



uk^Y u^k^VX^ ( \B\k 

■ arctan 



^\B\kJ 

Uk — uk"^ 
Uk + vk"^ 

arctan 



(|B|A:) 

' ujk + vk 
IBIA; 



3 

2\2 



Vv^A^ 



{uk + vk'^f^ujl + (A + v)k'^ujk + Az/A;4 
WW 



IBIA; 



^jul + (A + v)k'^ujk + \vk'^, 



(42) 



n,(A;) 



XeY X^k^Vi^Xk* ( \B\k 

arctan 



+ 



JB|A;; 

Xk - A/c^ 
Xk + AA;2 

arctan 



(|B|A:)2 



V^^aF 



'Xfc + A/tn ' {Xk + Afc2)yxi + (^ + A)A;2xfc + uXk^ 



IBIA; 



(|B|A;)= 



IBIA; 



'xl + {^ + ^)k\k + J^Xk\ 
With (24), (33), (38) and (39), we finally have an expression for a 



\ r+oo 

a = — / dk 
2 Jo 



^^ A,{k)U,{k) - ^^A,{k)U,{k) 



X + u 



Z/ + A 



■X 



(43) 



(44) 



4. The Dependence of a on Kinetic Reynolds Number R^, Magnetic Reynolds 

Number Rm and Prandtl Number Pr 

Before we discuss the dependence of a on different dimensionless parameters of 
the turbulence, we explain the physical meanings of the notations introduced in above 
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derivation. C{k) in relation (40) was first introduced in the correlation spectrum of zerotli 
order velocity field, i.e., equation (36). It is the helical part of the correlation function of 
velocity field v*^'^-'. In fact, one can easily obtain the following relation 

f 87rk*C{k)dk =< v(°) ■ V x v(°) >, (45) 

where < ■ > denotes spatial average. This shows that Siik'^C^k) can be regarded as the 
spectrum of kinetic helicity. In parallel, we may regard 8TTk'^H{k) as the spectrum of the 
current helicity of magnetic field according to the following relation 

f 8nk^H{k)dk =< b(°) ■ V X b(°) > . (46) 

l/cok is the de-correlation time of v*^°) at scale ~ 1/k, and 1/xk is the de-correlation time 
of b*^°) at scale ~ 1/k. These two parameters are statistical properties of MHD turbulence 
of v*^°) and b*^°) . They are independent of B but certainly affected by the external forcing 
term, f , as v*^°^ and b*^*^) are driven by this force. We believe that these two parameters 
represent the effect of turbulent nonlinear interactions between different modes of v*^°^ and 
b*-^-', such as the nonlinear cascade of energy from the forcing scale ~ /q down to dissipation 
scales. Note that in our model, ujk and Xk are functions of the magnitude of wave vector 
k, not its direction. This is because we made the assumption that the MHD turbulence 
of v*^°) and b*^*^) is isotropic. This is certainly a simplification. Because the turbulence is 
driven by external force f , if f is not isotropic, anisotropic structures of v*^*^) and b*^°^ will 
develop in turbulence, and the a-effect cannot be represented by a parameter as shown in 
(24). Rather, we must write a as a tensor. If external force is not strongly anisotropic, 
equation (44) may be used to approximate the isotropic component of the a-effect. For 
strongly anisotropic turbulence, the a calculated in above section only gives the a-dynamo 
effect from the isotropic parts of the velocity and magnetic field, which may play much less 
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iinportant role than the anisotropic parts.0 

If we consider the special case that molecular viscosity is the same as magnetic 
resistivity, A = z/, we have the following expression 

1 r+oo 

a = -- dk[T,,k87reC{k)U,{k) - n^kSnk^ H {k)Ilb{k)i (47) 

2 Jo 



2 Jo 
where 



1 1 , 

Tv,k = — ■,n,k = — , (4 J 

^k Xk 



For fixed, negligible but non-zero, values of z/ = A and finite values of r^ ^ and r^ y^, in the 



^Another concern of such isotropy assumption, is that recent work on MHD turbulence 
(see Moran 2000 and references therein) shows that anisotropic small-scale structure of v'^^-' 
and b^^-* develops within the inertial range; therefore the direct cascade of both kinetic and 
magnetic energy will not be isotropic. The numerical simulation by Moran(2000) shows that, 
if MHD turbulence at the forcing scale is isotropic, small anisotropic structure only develops 
at scales five or six times smaller than the forcing scale. That is, for the MHD turbulence 
of v*^*^) and b*^*^) within the scale range [e/o, e~^^o] ^oi e ~ 5, v^^^ and b^^^ are still isotropic to 
a large degree. But for scale / ^ e^^^o? because two initially adjacent magnetic field lines 
can stay correlated after they wander a distance of /, small-scale anisotropic structures will 
develop according to the anisotropic energy cascade mechanism (see Goldreich and Sridhar, 
1995). However, one must notice that such anisotropy is defined with respect to a magnetic 
field b^°)(/') for certain scale /' > / at certain instance. h^^\l') itself is isotropic when averaged 
over several turbulence outer scale Iq, or when averaged over many eddy turn-over time of the 
turbulence. Because the a effect considered in this work is an ensemble averaged quantity 
over the scale L ^ Iq, we believe the anisotropy of MHD turbulence at very small scales will 
not change our final results qualitatively. 
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limit oi B ^ 0, because 

oo /_ 
PiTC.tPiT\(f) = Y^ — 

2k + 1 



oo / 1 \kc2k+l 

arctan(0 = E o/ , i ' ^^ ^ ^ 0' (^9) 

fc=0 



we have 



n„(5) ^° ^. (50) 



With (38) and (39), we have the a effect in the hmit of small B as 

a^-^ (r, < v(°) ■ V X v(°) > -r, < b(°) ■ V x b^") >) . (51) 

In the limit of r^ = Tb, the same result has been obtained by Field et al. (1999). Here t^ 
and Tb are the typical de-correlation times of velocity field and magnetic field, respectively. 
Note that the current helicity term, |r5 < b*^*^) ■ V x b*^*^^ >, cannot be neglected. It appears 
naturally when v and b are treated on an equal footing, which was not considered by 
kinematic dynamo theory. 

We expect that a-effect is quenched if B is moderate or large, i.e., B > vq where Vq is 
the turbulent velocity at scale Iq. This can be seen from (47) under the conditions that: 
(a)-B > Vq, (b)z/(= A) is small, and (c)the de-correlation time is smaller than the eddy 
turn-over time but larger than the Alfven time ~ 1/kB. Because 

TT 

arctan(^) — > — , for ^ ^ oo, (52) 

with conditions (a) and (b), we have 

/veV vV^M* I \B\k \ fvk^V 

[■m) ~ ~mkr ""'"' [7m) ^[-m) "'' 

for small u, A and moderate or large B. Condition (c) requires that: (d) uJk/ Bk < 1. Thus 
with (d), (a) and (b), we have that up to the second order in uj^/Bk 



1^ |B|A; ) (|B|A;)3 ^ 

arctan ( , '^'^ \ -^ (^\ \ (54) 

V^cu2 + (A + z/)A;2a;fc + Az/A;V V|B|A;; 



a « — =2 / dk[T^,kS7rPC{k) - n,kSnk^H{k)], (56) 
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Also, to the second order in Xk/Bk, we have 

(Xk+j^y _ iXk + ^ky^xl + (A + iy)k'xk + A^ 
I |B|A: ) (|B|A;)3 ^ 

arctan ( ^ '^'^ ] _ ( ^] \ (55) 

With these approximations, we have that under conditions (a), (b) and (c), 

\ r+oo 

2B^ Jo 
where C{k) and H{k) are related to kinetic hehcity and current hehcity through (45) 

and (46). Note that this resuh is independent of z/(= A). (56) shows that a falls rapidly 

as B increases. To further relate a with the statistical properties of MHD turbulence, 

we introduce a vector u*^°) that is the vector potential of velocity field v*^*^) through 

V X u*^°^ = v*-"-*. Let a*^*^-* be the vector potential for h^^\ It is easy to find the relation 

/ 87reC{k)dk =< u(°) ■ v(°) >, (57) 

and 

/ 8nPH{k)dk =< a(°) ■ b^") > . (58) 

The right hand side of equation (58) is the magnetic helicity of the turbulent field h^^\ 
With (57) and (58), we may re- write (56) as 



/<u(o).v(o)> <aW-bW>\ , , 

a = - { =2 =2 • (59) 

V 2B\, 2B\, I 

(59) relates a-effect with the dynamics of magnetic helicity of turbulence. Because magnetic 
helicity is a conserved quantity in ideal MHD, a coefficient derived above will be constrained 
by the dynamics of magnetic helicity. 

Finally, we consider how non-unit Prandtl number affects a dynamo. To aid our 
discussion, we introduce the following quantities: 

RmAk) = ^, (60) 



ReAk) -- 


Z/P' 


Rmjk) : 


Xk 
AF' 


ReAk) -- 


Xk 
z/F' 


pm = 


B|A; 
^k 


Mk) = 


\B\k 
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(61) 
(62) 
(63) 

(64) 

^''^ 
These parameters are determined by the statistical properties at different scales of the MHD 

turbulence of v*^°) and b*^°\ They are all functions of k. Each of the first four quantities 

denotes the ratio of the kinetic or magnetic dissipation time to the de-correlation time of 

velocity field or magnetic field, therefore is similar to the definition of Reynolds numbers, 

which are the ratios of nonlinear term to the dissipation term in the momentum equation 

or the induction equation. f3i and (32 measure the strength of large-scale magnetic field in 

the units of v{k). Let t^ be the eddy turn-over time of the turbulence at the forcing scale 

/o- For these parameters at scale Iq = l/ko, with the experimental result (15), we have their 

relations to the conventional Reynolds numbers as 

RmAko) = Rm— ~ (3 ~ 5)Rm, (66) 



r, 



V 



ReAko) = Re- ~ (3 ~ 5)i?e, (67) 

RmAko) = Rm- ~ (3 ~ 5)i?„, (68) 

n 

ReAko) = Re- ~ (3 ~ 5)i?e, (69) 

n 

Piiko) = ^^ ^ (0.2 ~ 0.3)^, (70) 

P2{ko) = ^? ^ (0.2 ~ 0.3)ffl. (71) 

Vq tv Vq 
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Next, we express the general result of (44) using these statistical quantities of the 
MHD turbulence. We first write 

POO 

a = akdk, (72) 

Jo 

where 



ttfc 



^ -A,(A;)n,(A;) - .^^Ab{k)Uk{k) 



1 + Pr ^ ' ^ ' Pr + 1 
Here Pr = u/X is the magnetic Prandtl number. Define ZM^k) = ^_^p^ Ily{k) , 
ZN{k) = 3^Tl,{k)- then 



(73) 



ftfe = - [A^(A;)Zm - ^h{k)ZM] 



(74) 



With those quantities introduced in (60)-(65) and the definitions of Zm and Z^-, we have 



Zmik) 



+ 



1 + Pr 



Re,ui — 1 
Re,uj + 1 



arctan 



1 



ixp 



iPlRe,u 



arctan 



l/Pr 



lY {Re,^ + ifjRl^ + (1/Pr + l)R,^^ + 1/Pr 



PlRe,u 



iPlRe 



■X 



PlRe,u 



i?2^ + (l/Pr + l)i?e,.+ l/Pr^ 



(75) 



^iv(A;) 



+ 



2Pr 

Pr + 1 

J^m,x ' 

arctan 



,P2Rm.,^ 



arctan ■ 



{(^2Rm, 



'Pr 



R^,x + 1 V (^-.x + l)V^^.x + (^^ + l)^™-x + ^^ 



lj2Rm,x j 

P2Rm,x 



{l32R^,xY 



X 



P^^ + (Pr + l)P^,^ + Pr^ 



(76) 



If we compare (74) with the classical result given by (51), we find that the non-unit Prandtl 
number will modify the a coefficient through the complicated terms Zm and Zjy. For 
Pr -^ oo, we expect that Zm -^ 0; therefore, the contribution to a from the velocity field 
will be small. But large Pr may not switch off the contribution of magnetic field to a 
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because Z]^ can remain moderate for Pr -^ oo. On the contrary, for Pr -^ 0, Zjv -^ 0, 
so the main contribution to a comes from A^,. To further illustrate this, we plot Zm and 
Zn at different Pi and P2 in Figure 1. Zm is bounded above by | and decreases as l/Pf. 
When Rm,uj and j3i are fixed, Zm achieves its maximum when Pr falls in a certain interval. 
Outside the interval, Zm falls sharply to 0. Zj^ has similar properties (in fact, Z]^ and Zm 
are symmetric with respect to Pr = 1). In Figure 2 we show other views of Zm and Zn 
as functions of -Re,a; and Rm,x when f3i^2{= 0.0001) and Pr are fixed. The curve for Zm 
with Pr = 1 is already shown in Figure 1 of Field et al. (1999); however, as Pr increases 
or decreases from 1, the distribution of Zm as a function of i?e,w changes dramatically: 
for small Pr like 10~^, Zm stays at its maximum value | for all -Rg.t^ ^ 10^, compared to 
Zm — for all Re^uj k, 10^ when Pr = 10^, and Zm — | for all i?e,a; k, 10^ when Pr = I. 
But we must consider the contribution from Z^v, too. In panel (b) of Figure 2 we plot 
a few curves of Zjv as a function of Rm,x- Close investigation of panels (a) and (b) in 
Figure 2 shows that Zm and Zn may complement each other such that if the helicity terms 
|A^,| ~ |Ai|(see equation (74)), a can maintain a non-trivial value for a large range of Pr, 
and this is shown in Figure 3, where we plot the regimes where Zm > 0.1 or Zn > 0.1 for 
A = /52 = 0.0001. The vertically shaded area by solid lines is the regime of Re^w and Rm,Lu 
where the contribution to a from velocity field, Zm, is no less than 0.1 for f3i = 0.0001, 
while the horizontally shaded area by dash-dot lines is the regime of i?e,x and Rm,x where 
the contribution to a from magnetic field, Z^, is no less than 0.1 for (32 = 0.0001. The 
case with stronger large-scale magnetic field is shown in Figure 4 with /5i = /52 = 1. An 
{Re,uj, Rm,uj, -Re,x? Rm,x) Combination corresponds to two points on the plane of Figure 3 or 
Figure 4, one for {Re,ui, Rm,uj), the other for {Re,x, Rm,x)- ^^ either of these two points falls 
within its corresponding shaded area, or both do, we may obtain a non-vanishing a which 
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by (74) satisfies 

\a\ ;^ min{0.05|A^|,0.05|Aft|,0.05||A^| - IA5II} (77) 

at /5i = /52 = 0.0001. 

5. Discussion 

In our derivation for a, tlie de-correlation time for velocity field at scale 1/k, l/uv^k, 
can be different from the de-correlation time of magnetic field, l/xb,k, at the same scale. 
These two parameters are the statistical properties of the MHD turbulence, and can be 
measured in numerical simulations or experiments. Under different circumstances, these 
two parameters can be different. They are affected by molecular diffusion and magnetic 
diffusion, i,e., z/ and A, by the properties of external forcing term f, and by the nonlinear 
interactions between different modes of velocity field and magnetic field at different scales. 

We also generalize the derivation of the dynamo a-coefficient to include non-unit 
Prandtl number. Prandtl number is an important physical parameter. For many 
astrophysics systems, Pr 7^ 1. For example, in the solar convection zone, Pr ~ 10^® — 10~^ 
(Childress and Gilbert, 1995); while in the partially ionized warm gas of the interstellar 
medium Pr ~ 10^^ (Kinney et al. , 2000). Many past and recent discussions on a dynamo 
assume unit Prandtl number, which may not appropriate for real astrophysics systems. Our 
analysis of non-unit Prandtl number case shows that, for astrophysical systems such as solar 
convection zone where Pr can be very small, the contribution to a dynamo from kinetic 
helicity can be much stronger than the contribution from current helicity. In this limit, our 
result can be reduced to the classical result of kinematic mean- field theory, i.e., the current 
helicity term does not play important role in determining a. This can be seen from (73) by 
setting Pr -^ 0. However, for the opposite limit, i.e., Pr -^ 00, the contribution to a from 
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current helicity can be much more important than that from kinetic hehcity. This means 
for astrophysical systems such as warm interstellar medium, dynamo a effect will strongly 
depend on the current helicity of magnetic field of the system instead of the kinetic helicity. 

Another possible feature of a dynamo considered in this work is that, regardless of 
large or small Prandtl number, there is no simple relation as 

a oc ^ (78) 

1 + r^b'IvI 

for Rm ^ 1. Rather, nonlinear interactions of the turbulence velocity field v*^°) and b*^'') 

will lead to finite de-correlation times of v' and b' so that even for Rm,Re ^ Ij the 

electromotive force is not significantly reduced. This can be easily seen from the plot in 

Figures 3 and 4. a calculated in the shaded areas (horizontal or vertical) is not dramatically 

quenched. Such shaded areas span large parameter space of {Re,u},Rm,u},Re,xyRm,x)- Hence, 

we believe that for large kinetic and magnetic Reynolds numbers, the following relation is 

approximately true 

a oc =^— (79) 

B /v^ 

for moderate or large B (see expressions (56) or (59)). Here, for Pr ~ 1, 
ao = -y < V ■ V X V > +^ < b ■ V X b >; for Pr < 1, ao = -^ < v ■ V x v >; and for 
Pr ^1, Q;o = y<b-Vxb>. Relation (79) is close to relation (9), which Kraichnan 
proposed in 1979. Indeed, our model is closely related to, but not exactly the same as, 
Kraichnan's oversimplified model. A detailed comparison between our model and the model 
by Kraichnan has been given in appendix B. 

The derivation of dynamo a effect in this work should be considered as an extension 
of the kinematic MFE to the dynamic regime. We use two-scale approach, but we treat 
the velocity field and the magnetic field on equal footing. Our assumption of short time 
de-correlation (see relation (14)) has been proved by both experiments and numerical 
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simulations. There are, however, the circumstances that condition (14) is not vahd, such as 
the small-scale magnetic field amplification problem considered by Kulsrud and Anderson 
(1992). In their work, when magnetic field is weak so back reaction can be ignored, the 
de-correlation time of magnetic field can be long compared with the line stretching time; 
therefore, the nonlinear terms must be kept. Kulsrud and Anderson (1992) did exactly that. 
They assumed that the smallest time in the frame of their work is the turn-over time of the 
smallest eddies, Tmax, and the dominant effect is the line stretching by turbulent eddies. 
But their results can be valid only for small back reaction. In our model, back reaction can 
be strong, and the smallest time is the de-correlation time of the B-dependent v' and b'. 

It is also possible to incorporate into our model the nonlinear effects from those 
nonlinear terms that we omitted in section 2, if the de-correlation times of velocity and 
magnetic fields are long compared to eddy turn-over time. To aid our discussion, we focus 
on the case that B is strong. In this case, the MHD turbulence can be treated as a bath of 
Alfven waves. In appendix B, we argue that in the limit of strong large-scale magnetic field 
B, v' and b' can be regarded as {v*-"-*, b*-°-'}-driven, Alfvenic wave-like motions. Because 
v*^*^^ and b*^°^ can be de-correlated within Tdcor < Teddy, the {v*^°\ b*^°^}-driven v' and b' will 
also be de-correlated within Tdcor < ^eddy Our understanding of MHD turbulence is as 
follows: for statistically steady state of the MHD turbulence, the net effect of the nonlinear 
interactions between different modes of such wave-like motions, is that the kinetic energy 
and the magnetic energy will cascade from forcing scale down to dissipation scales and be 
converted to heat. Therefore, to incorporate the nonlinear interactions between different 
modes of Alfvenic wave-like quantities v' and b' into our model, we approximate the 
nonlinear terms with linear, effective turbulent mixing terms, and re- write equation (17) 
and (18) as 

(-lu + ue + andk') v'^ + indk% = ?(k ■ B)(6;. + Sf ), (80) 
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'-IUO + AA;2 + r]ndk'') h'^ + Wndk'^v'j = ?(k ■ B)({;; + vf). (81) 

Here we approximate the nonlinear terms in the (Fourier transformed) equations for v' and 
b' with (Jndk'^'^' + Ind.k'^^' and rjndk'^h' + Wndk"^^' (see Appendix B for more discussions). 
All of the nonlinear damping factors, and, Ind, Vnd and zund, are functions of k. These 
damping factors can be calculated with closure theory (Biskamp, 1994; see also Chen and 
Montgomery, 1987, or Kraichnan, 1979). In the limit of large B, compared with terms 
involving ih ■ B, the nonlinear damping factors '-fnd in (80) and Wnd in (81) have significant 
contributions to the damping of turbulence only near the plane that is perpendicular to B; 
therefore, their contributions to the integrations with respect to ^ = Z(k, B) are small and 
we omit them to get the first order approximation. Under such considerations, molecular 
viscosity u and magnetic diffusivity A in our calculations for a in section 4 must be replaced 
by z/ + and and A + rjnd, respectively, and the magnetic Prandtl number in expression (73) 
has to be replaced by 

^ + 0"nrf(^) / N 

and the ak in (73) must be re-written as 
1 



^k,nd „ 



^ A,(A;)n,,„,(A;) - „^^''"f ^\ A,(A:)n,,w(A:) 



l + Prnd{k) "' ' "'""' ' Prnd{k)^\ 
Here Hj, „(i(A;) and Hf, ^(^(A;) can be obtained by replacing z/ by z/ + and{k) and by replacing 
A by A + i]nd{k) in expressions (42) and (43). The approximation to nonlinear effects 
with and{k) and i]nd{k) are simple, straightforward, but crude. Detailed analysis of the 
exact forms of and{k) and rjndik) with closure theory or other analytic methods can be 
complicated and is beyond the scope of this work. 

It has been realized that the conservation of magnetic helicity is related to the a 
dynamo process. Many authors in the literature have suggested that for large magnetic 
Reynolds number, a effect is suppressed. Our model does not consider the conservation 
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of magnetic helicity. This is because we consider the a effect as a resuh of the nonhnear 
interaction between homogeneous, isotropic MHD turbulence driven by an external force 
and a large-scale magnetic field, B. In our model, the external forcing term has the freedom 
to drive turbulence that will interact with B to produce a non-zero a dynamo in large Rm 
limit, while in the work by Gruzinov and Diamond (1994, 1995, 1996; see also Zeldovich, 
1957), the velocity field is simply assumed to be given and no external energy sources are 
considered. This is not true for astrophysical systems where "...kinetic energy is constantly 
being injected hydrodynamically [at the forcing scales]." (Kulsrud, 1999). Note that the 
result that a effect should be quenched due to the conservation of magnetic helicity in the 
large Rm limit is obtained by assuming a steady state and a closed system. As Blackman 
and Field(2000) point out, the assumption that a system is closed forbids a non-zero net 
fiux of magnetic helicity to flow into or out of the system. In real astrophysical systems 
like the Sun, the boundary is open, and magnetic helicity can flow through it. Moreover, 
the assumption of steady state is not valid in many astrophysical systems where transient, 
impulsive activities, such as solar flares, can happen(Chou and Field, 2000). With no 
assumptions of closeness or stationarity of the system, we believe the a dynamo can 
operate under the nonlinear interaction between B and the MHD turbulence, thus can be 
determined by B and the statistical properties of the turbulence. A more thorough, though 
complicated, model of a dynamo than ours should incorporate the boundary effects into 
the discussions. 

6. Conclusion 

We generalize the derivation of the dynamo a-coefficient to include non-unit Prandtl 
number and different de-correlation times of velocity/magnetic field. Our formula gives 



-28- 

q; as a functional of the statistical properties of b'^*'^ and v'^*'\ and it is also a function of 
several parameters introduced in section 4. We confirm the fundamental results of Field, 
Blackman and Chou(1999), that a does not vanish at very high Reynolds numbers, but 
rather obtains a finite value. There are two parts in the expression of a: one is from 
velocity field, the other is from magnetic field. Non-trivial contribution of either of these 
two parts can be achieved only when the Prandtl number Pr falls in certain regimes, but 
the total contribution to a can be insensitive to Pr if kinetic helicity and current helicity 
are comparable. 

I am grateful to Prof. George Field for many stimulating discussions. I also benefit 
from my discussions with E.G. Blackman, A. Brandenberg, B. Chandran, S. Gowley, R. 
Kulsrud, D. Balsra. 

A. Calculation of Electromotive Force 

We start from the calculation of correlation function of velocity field and magnetic 
field. If we plug (19) and (20) into (25) in section 3, we have 



cikdcue^^'^-''-'^*) (t;Kk', cj')&i(k - k', cj - cj 



(t;/(x,t)6j(x,t)) = / / d\lduj 

d\^dJ{\\ - C(-k', -J) - C(k',cj') + 



c(k',c.')c(-k',-c.')]p,';°(-k',-c.') + 

[a(-k', -J) - c(k',c.oa(-k', -c.')]Qi?(-k', -j)\ + 

[e.(k', J) - C(-k', -c.')C.(k', c.')]^g°(-k', -J)\ + 

[e.(k',c.')a(-k', -j)\p'^iy,j)}. (Ai) 
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Note that [1 — C(~k, —u) — ((k, u) + ^(k, uj)({k, —uj)] and [^jy(k, uj)^x{—^, ~^)] are even 
functions of k, while [^a(— k, — c<j)— ^(k, cij)^a(— k, — c<j)] and [^^{k,uj) — ({—'k,—uj)^^(k,uj)] are 
odd. Because of the homogeneity and isotropy of the turbulence, (5['j'(— k, cu) = Q^^(k,uj); 
therefore the antisymmetric tensor (5^°(k, cj) — Q^^{\i.,uj) is an odd function of k, as are 
P;^°(k, u;) — Pj'i^(k,u) and Sij^(k,uj) — S'°f(k, cj), so many terms in (Al) vanish. By these, 
we have the electromotive force (v x b)^ = M„ + Nn- Mn and Nn are given in equations 
(34) and (35) of main text. 

With the definitions of ^x, ^u and ( in (21), (22) and (23), we can rewrite M„ and Nn as 
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(A4) 



(A5) 



In the following, we assume the square roots in (A4) and (A5) are real (otherwise, see Chou 
and Fish, 1999). Under the homogeneity and isotropy assumptions made in section 3, M„ 
can be computed as 



M^='-^l rfkX°°(k) 



(k Bl 



dipT{ip) 



{uJ — uji){uj — uJ2){u; — ujs){uj — UJ4) 

(A6) 



Contour integration of (A6) gives 

1 



M„ = — I rfke„,,XgO(k)[-(k ■ B)]{2m)[U - V] 



(A7) 
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where 
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Because Yl, Y2, Wl, W2 are not functions of ip, U and V can be calculated as 
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where 



7 



(cc^i - txe){ujj + u^k^) {uj2 - ixe){uji + u^k^) 



Because of the following five identities: 
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UJi — UJ2 = UJs — UJ4 

with (A4) and (A5), we have then 

Xu^k"^ LUkX 



-(A-z/)2A;4 + 4(k-B)2, 



(A18) 

(A19) 
(A20) 
(A21) 
(A22) 



(A + v)ui)k X + V 
Xu^k^ : 



ujI + {X + u)k'^uJk + Xuk^) +(A;|B|cos 



{X + iy)ujk Xuk^ + {k\B\ cosOy' 

Here 6 is the angle between vectors k and B. With (36) and (37) in section 
3 of main text, we have that e„ijX['P(k) = ienij€ijkkkC{k) = i2knC{k), and 



(A23) 



(^nijO'^j'iX) = i^nij(^ijkkkH{k) = i2knH{k). Let p = ^ul + {X + u)k'^iJk + Xuk'^, 
y = A;|B| cos^, and define Q{y) as 



e(y) 



Xu^k^ - A, 



UJt 



Xv'k 



2lA 



(A + U)ujk p2 + y2 (^\j^ ^^^^ \j^j^2 j^ y2 



(A24) 



Then M„ can be rewritten as 



M„ = -47r/ dkk'^C{k)Bn cos"^ 9 d{cos9) Q{k\B\cos9). (A25) 
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Note that 



cos2^rf(cos^)0(A;|B|cos^) 



\uol - Wk^ 1 



(A + z/H {k\B\fJ-mp^ + y 



fe|B| y2 



,2-^^ 



+ 



Xv'k 



2lA 



fclBI 



y 



-dy 



(A + v)uJk {k\B\f J-k\B\ {y/X^y + y2 
__^n,(A;) (A26) 

UJk A + U 



where 



Ilv{k) 



+ 



uk 

Mk) 

Uk — i^k"^ 
ujk + vk'^ 

arctan 



2 \ 2 , ,2 l4 



u^k'^VXuk^ ( \B\k 

— ,, — , ^^ — arctan , 
(|B|A;)3 \^hJB 



'oJk + iyeV i^k + i^k'^fJ^k + (A + y)k^ujk + Az/A;4 



IBIA; 



(|B|A;)= 



■X 



IBIA; 



Let 



cj| + (A + u)k'^ujk + \uk^ 
^nk^Cik) 



therefore we have 



In parallel, we have 



where 



A„(A;) = 
Ab{k) = 



\ r+oo 

M„ = -- / dk 
2 JO 



^k 
Mk^H{k) 



Xk 



2A 



X /■+00 

Nn = - dk 

2 Jo 



2u 



A,{k)U,{k) 



Br 



Z/ + A 



A,{k)U,{k) 



Br 



(A27) 

(A28) 
(A29) 

(A30) 

(A31) 



n,(A;) 



+ 



"/ 


xk^y 


v 


Bk) 


Xk 


-\k^ 


Xk 


+ AP 



A2A;4v/^;aF / |B|A; 

arctan 



{\B\ky 



ViAk^, 



Xk + Xk'Y iXk + XkyJxl + (^ + \)k'Xk + y\k' 



arctan 



IBIA; 



mkf 



■X 



IBIA; 



xl + {u + \)k^Xk + u\k\ 



(A32) 
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We finally achieve an expression for the electromotive force 

B. (A33) 



\ r+oo 

< V X b >= — / dk 

2 Jo 



^^ -A.(fc)n.(fc) - -^A,{k)U,ik) 



X + u "' ' "' ' Z/ + A 



B. Alfven wave-like properties of {v', b'} in the presence of strong B 

In the presence of a strong large-scale magnetic field B, MHD turbulence takes the 
form of random Alfven waves propagating along B (Kraichnan 1965). This can be seen 
from (12) and (13). Indeed, if we re-write (12) and (13) as 

9iv' = B- Vb' + {z/VV + B- Vb(°)}, (Bl) 

aib' = B- Vv' + {AVV + B- Vv^}, (B2) 

we find that without the terms in brackets, equations (Bl) and (B2) have solutions of linear 
Alfven waves. These Alfven waves are damped through the dissipation terms. External 
energy sources, represented by the forcing term f in the momentum equation (11), will drive 
a MHD turbulence represented by {v^'^\ b'^'^^}, which in turn stretch the large-scale magnetic 
field B to drive the Alfven-wave-like patterns of {v',b'}. Such physical processes and the 
damping due to {z/. A} are represented by the terms in bracket in equations (Bl) and (B2). 
Therefore, in the presence of strong B, {v', b'} are Alfvenic wave-like motions that are 
driven by {v^^^ b^^^} and damped by {z/. A} terms. To the degree that the condition (14) is 
not exactly obeyed, such {v^°\ b^^^j-driven, {z/, A}-damped, Alfvenic wave-like motions can 
be further damped by nonlinear wave- wave interactions, as discussed by Kraichnan (1979). 
In the limit of large 5, if we use linear terms Ondk^ ■, Indh^ ■, Vridk"^ and Wndk"^ to approximate 
nonlinear damping term^, we may write, in Fourier space, the following two equations for 



^There terms are, using our notations in section 2 of main text, the following: b' ■ Vb', 

b' • Vb(o), b' ■ Vv', b' • Vv(o), v' • Vb', V ■ Vh^^\ V ■ Vv', v' ■ Vv^o), b^ . Vb', b^o) . Vv', 
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v' and b' 

{-zu + uk^ + a„,A;2) v'^ + ^^,k%'^ = ^(k ■ B)(b'^ + bf), (B3) 

(_,^ + Xk' + r^„,A;2) &;. + w^,ev'^ = ^(k • B){v'^ + {;f )). (B4) 

These two equations are modifications to equations (Bl) and (B2). Following Kraiclinan 
(1979), we consider single mode solutions of (B3) and (B4): 

b'(k) oc (k • b) e-'^^^'-'^-^e-^'^''^'''' (e''^''-'=*b(°)(k) - e'<^-'=V(°)(k)) , (B5) 

v'(k) oc ± (k ■ b) eT^^-^'=*Tir„g-r„(fc)fc2i (^^iK,,%m ^^^-^ _ e^<^^.^V(°)(k)) . (B6) 

Here Va = -B is the Alfven speed in the unit of {AirpY''^. r^ and r^ are phases determined by 
initial conditions. r^y){k) are linear combinations of the damping effects X-iV^Undiilndilnd 
and Wndi i-e., damping effects from both nonlinear mode-mode interactions and physical 
viscosity/diffusivity. (j)\j\, = gJk + i^«,A; is the phase of v*^°\ while (Pik ~ 9hk + iX6,fc is the 
phase of b'^°\ Note that depending on the initial condition, gv{b),k may not be equal to 
9v(b),k- Relations (B5) and (B6) tell us the following points: (1) v' and b' are Alfven waves 
in character; yet, (2) these Alfven waves are driven by {v'^'^\ b*^°^}, which are in turn driven 
by external forcing term f; (3) these {v*^°\ b*^°^}-driven Alfven waves are damped in three 
ways: (A) molecular viscosity u and magnetic diffusivity A, and (B) nonlinear interactions 
between different modes of {v', b'} approximated by parameters (TndjVndjInd and 'CCnd, 
and (C) the de-correlation of {v*^°\b*^°)} approximated by cu^^k and Xb,k- The a dynamo 
calculated from (B3) and (B4) in the limit of large B is 



oo 



a = ak,nddk (B7) 







where ak,nd is given by (83). 



v(°) ■ Vv', and v(°) ■ Vb'. 
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Kraichnan (1979) demonstrated the a effect in the strong-field hmit where the 
turbulence consists of Alfven waves propagating along B. He considered the following 
linearized single mode helical standing wave, 

V = Mo cos(VaH) [sin(A;z), co^ikz), 0] , (B8) 

b = mq iwiiyA^t) [cos(A;z), — sin(A;z), 0] . (B9) 

He then calculated an a effect from such helical waves and found 

Here, I/7 is the typical nonlinear damping time.[| In the limit of large B, (BIO) can be 

approximated as 

a{k) ~ — =2~ (<u-v>-<a-b>), (Bll) 

2B T 

where V x u = v and V x a = b. r = I/7 is the de-correlation time. This result is the 
same as (59) of our calculation in large B limit, except the following three aspects: first, in 
our model the de-correlation times for v and b are different; second, the kinetic helicity 
and current helicity in our model come from v*^°^ and b*^*^); third, the nonlinear damping 
effect in our model comes from not only the nonlinear interactions between modes of Alfven 
waves as proposed by Kraichnan, but also the nonlinear interactions between {v*^°\ b^"-'}. 
The standing helical waves in (B8) and (B9) are solution of linear MHD equations with 



'^Note that 7 must come from the nonlinear interactions between Alfven wave modes of 
the type (B8) and (B9) at different wave lengths. This is because, given the single mode 
standing wave pattern like (B8) and (B9), the nonlinear terms in momentum equations, 
V X (V X v) and (V x b) x b, are both zero, so that velocity field can be damped only 
through molecular viscosity z/, and the magnetic field can only through magnetic viscosity, 
A. 
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special initial conditions. In our model, kinetic helicity and current helicity are determined 
by {v*^°\ b*^*^^}, which are driven by external helical forces. In many rotating astrophysical 
systems, such external helical forces can be buoyancy forces, or Coriolis forces. We believe 
the helical properties of the linear Alfven wave of Kraichnan (1979), i.e., relations (B8) 
and (B9), are determined by the MHD turbulence {v^°\ b^°^}, which is in turn driven by 
external helical forces (see also Moffat, 1978). Such external energy sources to the MHD 
turbulence is of vast importance to the dynamo action (Kulsrud, 1999), but have been 
overlooked in previous calculations of a dynamo. Our work is certainly an extension toward 
that direction. 
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Table 1. Notations and Their Physical Meanings 

Notation Physical Meaning 

V, A molecular viscosity and magnetic diffusivity 

Re,Rm kinetic and magnetic Reynolds numbers 

Pr — V jX magnetic Prandtl number 

a. a coefficient of dynamo theory 

B large-scale magnetic field, assumed constant 
v(o) ^ b*-*^) turbulent velocity and magnetic field that are independent of B 

v' , b' turbulent velocity and magnetic field that depend on B 

^0 = 1/fco size of energy containing eddies of the turbulence 

k magnitude of wave vector |k| 

ujkixk) de-correlation frequency of velocity field v*^°)(b*^°)) at k 

Ty{Tb) de-correlation time of turbulent velocity field v^"'(b''°^) 

Re,u,{k) = LJk/i^k^,^ (0.2 - 0.3)Re at fco 

Re^^ik) = Xfc/^^fc^ - (0.2 - 0.3)i?e at fco 

R„,.Uk) - ujk/XP,^ (0.2 - 0.3)i?,„ at fco 

Rm,xik) ^ Xk/Xk^, - (0.2 - 0.3)i?„ at fco 

Piik) = |B|/(c.fe/fc) 

Hk) = m/iXk/k) 

A^, = / A^{k)dk (XTy < v(o) • V X v(o) > 

Afe = / Ab{k)dk oc T6 < b(o) • V X b(") > 
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Fig. 1. — This plot shows Zm and Zn, defined in (75) and (76), as functions of Prandtl 
number, kinetic Reynolds number and magnetic Reynolds number at various values of (3i 
and P2 defined in (70) and (71). 



-42- 



(a)|3 =0.0001 




Pr=0.25 



SOOOOOOOOOOOOOijBOOOOOOOOOOOOOOOOOOOOOOOO Pr= 1 

-10 

Pr=5 



-0.1 



10 



10 



10 



10 



15 



10 



20 



e.co 



(b) P2=0.0001 



fl I I I I I I I I I I I I I I I I 




Pr=5 



OOOOOOOOOOOOO0OOOOOOOOOOOOOOOOOOOOOOOO Pr= 1 



Pr=0.25- 



-0.1 



10" 



10-^ 



10' 



10' 



10' 



m,x 



Fig. 2. — Zm{^ and Z^ih) at various values of Prandtl number. Here f3i = (32 = 0.0001. 
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Fig. 3. — The regimes where Zm > 0.1 (shaded by sohd hnes) or Zjv > 0.1 (shaded by dash-dot 
hnes) at various values of kinetic and magnetic Reynolds numbers. Here f3i = (32 = 0.0001. 



-44- 



20' 



18 



16 



14 



12 



110 



cc 



^ 8 



10 e,y, 



8 10 12 14 16 18 2 



"T r 



ii-2=-€-^=M 



"T 1 1 r 



(i 



::l 



i\ir 



;>;= 



.1 



^c 



18 



16 



14 



12 



10 ^ 

DC 

o 

8 o 



_2 I I I I I I I I I I I I _2 

-2 2 4 6 8 10 12 14 16 18 20 

|°gioRe 



e.Q) 



Fig. 4. — The regimes where Zm > 0.1(shaded by sohd hnes) or Zn > 0.1(shaded by dash- 
dot hnes) at various values of kinetic and magnetic Reynolds numbers. Here f3i = j32 = 1- 
Notice that the area covered by the shade is smaller than the shaded area shown in Figure 
3. This is due to the quenching of a-effect by large value of B. 



